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This paper discusses the aerodynamic behaviors of a gas mask canister with a complex inner structure 
and two porous materials in the filter layer and the activated carbon layer. The effects of the distribution 
and area of holes in the main sieve diaphragm and the thickness of the activated carbon layer on the 
pressure drop and the flow structure were determined using computational fluid dynamic (CFD) tools. 
The momentum loss of porous flow calculated by Forchheimer’s equation was added to the source term 
in the momentum equation. Streakline flow visualization was employed to observe gas flow structures 
within the empty canister and to identify the shortcomings of the prototype canister. Simulation results 
for the estimated inertial and viscosity parameters in Forchheimer’s equation agree closely with experi¬ 
mental values. The porosity of the canister for intake flows of 15-135 L/min causes the flow behavior to 
transition gradually from linear (viscous effect) to slight non-linear behavior (slight inertia effect). This 
study uses air age as an index of the time that air resides within the canister to displace the adsorption 
time of toxic gas. This approach conveniently elucidates overall filter capacity and the positions of dead 
zones in the activated carbon layer. The simulation results reveal that the channel design of the main 
sieve diaphragm dominates the aerodynamic behavior of the fluid within the activated carbon layer. 
Better hole distribution and a larger hole area correspond to a lower pressure drop, a smaller dead zone, 
and a higher adsorption time. The results in this study provide a valuable reference for designing channels 
in the main sieve diaphragm, and will be helpful in designing gas mask canisters. 

© 2009 Elsevier Ltd. All rights reserved. 


1. Introduction 

After the September 11 terrorist attacks on the World Trade Cen¬ 
ter in New York, a series of anthrax powder threats caused global 
anxiety. This type of threat caused an unexpected surge in gas mask 
sales. Improving the performance of gas masks therefore not only 
enhances personal safety, but also represents a huge commercial 
opportunity. The toxin-filtering function of a gas mask depends on 
the canister attached, which filters aerosols and absorbs toxic gases. 
Gas mask canisters typically consist of a filter layer and an acti¬ 
vated porous carbon layer that can efficiently block particles (e.g. 
anthrax), chemical agents (e.g. nerve gas), and germs (e.g. smallpox). 
The American National Institute of Occupational Safety and Health 
(NIOSH) and the US Army standard for general canisters both require 
that solid particles larger than 0.3 pm be filtered. Canisters are cate¬ 
gorized according to the efficiency with which they remove aerosol 
spray (toxic gas or droplet): 95%, 99% or 99.97% (Fischetti, 2002). Re¬ 
quirements for the next-generation M50 JSGPM mask developed by 
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the US Army state that weight and respiratory resistance must be 
50% lower than for the M40 series masks currently in service. 

This study analyzes a canister that has been marketed since 1970 
that has a filter capacity that meets the NIOSH 99.97% standard. 
However, its shortcomings are too heavy (330 g) and a high respi¬ 
ratory drag factor. At an inspiration flow rate of 30 L/min, the mea¬ 
sured pressure drop is about 200 Pa. The canister is 20% heavier than 
a similar product (Scott NBC22), and its pressure drop is 40% higher. 

The canister’s internal configuration (Fig. 1(a)) is very complex. 
The bracket helps support the filter layer, which comprises multi- 
pleated filter papers that trap large suspended aerosol particles. The 
main sieve diaphragm separates the filter layer from the activated 
carbon layer, and its structural strength alleviates the strain on the 
activated carbon layer during high stress movements. This design 
prevents gaps that can result in the loss of activated carbon parti¬ 
cles due to deformation, which in turn can cause canister failure. 
The main sieve diaphragm has several holes of various diameters, 
ranging from 4 to 7 mm. These holes allow gas to pass through 
the filter layer into the activated carbon layer. In the original de¬ 
sign, the outermost-ring holes in the main sieve diaphragm are 
sealed to prevent glue from migrating into the activated carbon 
layer during packing, as Fig. 1(b) shows. The type of glue used is 
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Fig. 1. (a) Dimensions of gas mask canister and components and (b) configuration 
of main sieve diaphragm. 


Table 1 

Main physical characteristics of the filter layer and activated carbon layer. 


Filter layer 

Pleated number 

36 

Material 

Polypropylene melt-blown non-woven 

Activated carbon layer 

Carbonization temperature (°C) 

1000 

Activation process 

H 2 0 

Particle diameter (mm) 

0.70-0.90 

Density (gm -3 ) 

0.40-0.50 

Porous volume (cm 3 g _1 ) 

>0.5 

BET surface area (m 2 g _1 ) 

> 1000 

Water content (%) 

<5 

Ash content (%) 

< 10 

Microporous volume (cm 3 g _1 ) 

>0.32 


polyurethane (PU); its function is to seal around the canister so that 
the activated carbon particles will not permeate into the filter layer. 
The activated carbon layer comprises activated carbon particles that 
absorb and filter toxic gases. The sub-sieve diaphragm prevents ac¬ 
tivated carbon particles from entering the canister outlet. The main 
physical characteristics of the filter paper and activated carbon par¬ 
ticles are given in Table 1. 

In common gas mask canisters, the mean inhalation flow rate 
varies from 15 L/min at rest to 135 L/min for intense exercise (Wang 


and Wang, 2003). Standard tests measure the pressure drop at a 
flow rate of 30 L/min. For the canister investigated in the present 
study, the entry Reynolds numbers (approx. 900-8100) calculated 
from the inhalation flow rates indicate that turbulent effects in the 
inflow should be analyzed. 

This study uses fluid dynamics to analyze the aerodynamic be¬ 
havior of the filter and activated carbon layers in the canister. These 
layers are porous media that cause a considerable pressure drop. In 
1856, Darcy first derived an equation to describe the momentum dis¬ 
sipation of a fluid flowing through a porous medium (Dullien, 1979): 


where AP is the pressure drop of the porous medium zone; L is the 
length of the flow direction; k is the permeability; p is the fluid vis¬ 
cosity, and denotes the superficial velocity of the fluid entering the 
porous medium zone. Darcy’s equation considers only the viscosity 
effect. A corrective convection term must be added to consider the 
inertia effect. This addition yields Forchheimer’s equation (Dullien, 
1979) 


A p 
L 


= a juV s + fipVj, 


( 2 ) 


where a is the reciprocal permeability of the porous material, 
or viscosity parameter. /? is usually called the inertial parameter. 
Forchheimer’s equation describes the relationship between the 
global pressure drop in the porous media zone and the superficial 
velocity on a macroscopic perspective. A higher superficial velocity 
corresponds to a greater pressure drop. The factors a, /?, L and V s 
determine the pressure drop. When the superficial velocity V s is 
very low or the Reynolds number is small, the inertia effect is far 
smaller than the viscosity effect. As the Reynolds number increases, 
the flow characteristics of the porous media gradually change from 
linear (Darcy’s equation) to non-linear (Forchheimer’s equation). 
This behavior differs from the general flow behavior, which involves 
an instantaneous change from laminar to turbulent flow (Dullien, 
1979). The non-dimensional form of Eq. (2) is 


'=s* + 1 - 


(3) 


where / = -Ap/LfipVj, and Re' = [3pV s /ocp. This equation is called 
the friction factor-Reynolds number form equation (Andrada et al., 
1999). This equation can be used to analyze linear to non-linear 
changes in porous media flow characteristics under various flow 
conditions. 

Aerodynamic behaviors, including pressure drops and flow struc¬ 
ture, are important considerations in the design of porous filters 
(Subrenat et al., 2000, 2003). Flow structures affect gas treatment 
processes within the filter by preventing preferential flow (which 
leads to absorbent failure due to overuse) and the formation of a 
dead zone (which leads to absorbent waste due to fluid not passage). 
Since the absorbent is opaque, the inner flow structures cannot be 
visually examined in terms of the size and location of preferential 
flow and dead zones. However, exploring the aerodynamic charac¬ 
teristics of a flow field using computational fluid dynamic (CFD) tools 
is helpful in analyzing the flow variables in porous media within a 
filter (Chen et al., 1995; Liu and Wang, 1996). The gas residence time 
in the filter must also be considered to ensure that the filter’s ab¬ 
sorption capacity is fully exploited. This study uses an air age index 
to determine whether or not overall filtering capacity is maintained 
after the canister body’s structure has been changed. A greater air 
age value corresponds to a longer air residence time. A lower air 
age value represents fresher air (Li et al., 2003). Air residence time 
distribution can therefore be used to optimize the geometry of in¬ 
dustrial filters (Baleo and Le Cloirec, 2000a, 2000b). This approach 
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conveniently elucidates overall filter capacity and the positions of 
dead zones in the activated carbon layer. 

Previous authors have used CFD tools to analyze porous filters, but 
their primary focus was to analyze the fluids flow through different 
arrays of the spatial microstructure of a filter media in 2- or 3-D 
simulations (Dhaniyala and Liu, 1999; Liu and Wang, 1996; Thomas 
et al., 1999). The geometrical design of particular filters, such as 
cylindrical pleated filters, was another purpose for using the pleating 
effect to design a new industrial filter (Baleo et al. 2000; Subrenat et 
al. 2000, 2001, 2003). The point which most previous studies have 
in common is that they focused on a single type of porous media. 
However, very little existing literature analyzes the aerodynamic 
characteristics of a gas mask canister. Therefore, the topic of the 
aerodynamic behavior of a gas mask canister with two kinds porous 
materials is worthy of further research. 

The purpose of this study is to evaluate the aerodynamic behav¬ 
iors of a gas mask canister in an attempt to decrease its pressure 
drag and overall weight. A streakline flow visualization reveals the 
flow phenomena of gas entering the empty canister as the outer¬ 
most holes on the main sieve diaphragm open and close. Observa¬ 
tions of this behavior were then compared to computational results 
to identify any problems with the prototype canister design. The 
pressure drops for the filter and activated carbon layers were then 
experimentally measured at different inlet Reynolds numbers. The 
Forchheimer coefficients in Eq. (2) are determined from measured 
values of pressure drop and velocity in either the filter layer or the 
activated carbon layer under different flow rates. The results for the 
porous media and compared to computational results to verify their 
accuracy and the effectiveness of the computational program. The 
aerodynamic behaviors of the canister were investigated for a range 
of inhalation flow rates by simulating different hole distributions and 
areas on the main sieve diaphragm, the thickness of the activated 
carbon layer, and the air age index. 

2. Problem 



Fig. 2. Mesh system of prototype canister (model Ai). 


where z is the air age means the local average residence time; Meff/ a r 
is the local actual diffusivity in turbulent model; p e jj = p { + p t (where 
\x e ff is effective viscosity, p/ is molecular viscosity and p t is turbulent 
viscosity) and cr T = 1 is the turbulent Schmidt number. 


2.1. Governing equations 


The governing equations herein include the continuity and mo¬ 
mentum equations, both of which obey the conservation principle 
(Subrenat et al., 2003), 

flm) = o. (4) 

CX[ 


d . . dp dzu d , —t—r. _ 

-{pu iUj ) = - — + ^ + m - + s,-. 

C/A j C/A j C/A j 


(5) 


where p is the gas density; is the velocity component in direc¬ 
tion i. The momentum Eq. (5) are steady, 3-D Reynolds-averaged 
Navier-Stoke’s equations, p is the pressure; zy is the viscosity shear 
stress tensor; g z is the acceleration due to gravity in direction i. The 
direction of the gravity term is perpendicular with the canister in¬ 
let, which is the axial direction (the opposite direction of z). When 
the flow is not laminar, the Reynolds stress term pu-uis related 
to the mean flow by the Boussinesq hypothesis. The source term 
Si describes the pressure gradient in the porous medium, and is 
given by Eq. (2) and is assumed to be isotropic. The inlet Reynolds 
number ranges from 900 to 8100, so this study adopts a modified 
low-Reynolds number k-s turbulent model (Jones and Launder, 
1973). The working gas is air, for which p = 1.225 kg nr 3 and 
p = 1.7894x1 0- 5 kg nrV 1 . 

The air age equation (Li et al., 2003) is 


2.2. Boundary conditions 

The canister inlet maintains a steady volumetric flow rate in all 
simulations. Since the inlet area is fixed, the calculated mean ve¬ 
locity is specified by the inlet boundary. A flow rate of 30L/min, 
which is regarded as the mean inhalation flow rate of normal adults 
in light motion, serves as the reference value for subsequent sim¬ 
ulations. The reference pressure is fixed at the inlet, with a value 
of 101,325 Pa. The outlet boundary condition is the pressure outlet 
boundary. The no-slip condition is assumed for all solid walls. For 
computational efficiency, only one-quarter of the full 3D model is 
constructed. 

2.3. Numerical method 

This study uses the FLUENT 6.2.16 flow solver, which is based on 
the finite volume method. The integral form of the governing equa¬ 
tions is discretized using an unstructured tetrahedral grid. The con¬ 
vection term is discretized with the first-order upwind scheme and 
the viscosity term is discretized with a second-order central differ¬ 
ential scheme. The SIMPLE algorithm (Patankar, 1980) is the solution 
algorithm for pressure-velocity coupling, and related discretization 
algebraic equations are solved using the TDMA method. 


2.4. Grid configuration and research matrix 


_d_ 

dxi 


(pm) 


_d_ 

dxi 


( fejf dl \ 

\ ^ fai) 


+ P. 


( 6 ) 


Since the filter layer material is very light and causes a much 
lower pressure drop than the activated carbon layer, this study 
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Table 2 

Research matrix. 


Model 

Parameter 




Hole distribution: outmost holes of main sieve 
diaphragm opened: O and closed: X 

Ratio of activated carbon layer thickness 
(note 1 ) 

Hole area effect: magnification of holes 
size (note 2) 

A 

Ai 

X 

1 

1 

a 2 

O 

1 

1 

B 

Bi 

X 

0.5 

1 

b 2 

O 

0.5 

1 

C 

Ci 

X 

0.25 

1 

c 2 

O 

0.25 

1 

D 

Di 

X 

1 

1.5 

d 2 

O 

1 

1.5 

E 

Ei 

X 

0.5 

1.5 

e 2 

O 

0.5 

1.5 

F 

f 2 

X 

0.25 

1.5 

f 2 

O 

0.25 

1.5 

G 

Gi 

X 

1 

0.8 

g 2 

O 

1 

0.8 

H 

Hi 

X 

1 

1.2 

h 2 

O 

1 

1.2 


Note 1: the ratio of activated carbon layer thickness for model A is 1. The activated carbon layer thickness of model B and model C is 0.5 and 0.25 of model A, respectively. 


Note 2: the hole size distributed on the main sieve diaphragm has a ratio of 1 in model A. The holes size for models D, G, H are 1.5, 0.8 and 1.2 times of model A, respectively. 



Flow Meter Terminal Box Mask 


Fig. 3. Schematic diagram of streakline flow visualization experiment. 


focuses on the activated carbon layer. To simplify the computa¬ 
tional domain, this study omits the effect of multi-pleated filter 
paper and treats the filter layer as a single layer. The medium of 
the two porous layers is assumed to be homogeneous. An unstruc¬ 
tured tetrahedral mesh with a complex structure and self-adaptive 
characteristics is adopted. Local grids for the inlet, outlet, sieve 
diaphragm holes, and the activated carbon layer are refined to 
accurately capture the flow patterns. The domain is divided into 
23 sub-domains, yielding a grid number of approximately 830,000 
grids for the prototype. Fig. 2 displays the domain mesh. Grid 
distributions for the subsequent models were generated based on 
the prototype mesh. All the results presented in this study are grid 
insensitive. The residual convergence threshold in all cases is less 
than lxlO -4 . 

The three most important parameters in determining the aerody¬ 
namic behavior of a gas mask canister are the distribution and area 
of holes on the main sieve diaphragm, and the thickness of the ac¬ 
tivated carbon layer. Table 2 presents the research matrix for more 
details. 


3. Results and discussion 

3.1. Flow observation and empty canister simulation 

Before adopting CFD simulations, streakline visualizations were 
used to observe the flow structures of an empty canister without 
filter paper and activated carbon particles. The inlet flow rate was 
30L/min. The experimental set-up included an argon laser, a smoke 
generator, a flow meter, and a camcorder, as Fig. 3 shows. The results 
reveal that a large backflow zone forms when the outermost holes 
are sealed. As a result, the fluid does not flow to the outlet easily, as 
Fig. 4(a) illustrates. Computational results reveal a similar outcome 
(Fig. 4(b)). Both the experimental and the simulated results demon¬ 
strate that the backflow zone disappears when the outermost holes 
are open and the fluid can flow out smoothly, as Fig. 5(a) and (b) 
shows. 

These results demonstrate that sealing the outermost holes in the 
main sieve diaphragm is a design flaw. When fluid flows actually 
through the activated carbon layer, the porous flow will not cause 
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Fig. 4. (a) Experimental image of activated carbon region (without activated carbon 
particles) at a flow rate of 30L/min (closed outermost holes on the main sieve 
diaphragm) and (b) simulated velocity vector of activated carbon region (without 
activated carbon particles) at a flow rate of 30L/min (closed outermost holes on 
the main sieve diaphragm). 


backflow in the peripheral zone. However, if the fluid cannot flow 
into this zone easily, a dead zone forms. In turn, this reduces the 
flow-through area and increases both the superficial velocity and the 
pressure drop. Computational simulations verify this inference. 

3.2. Porous medium parameters 

The filter and the activated carbon layers in the porous media 
are responsible for the majority of the canister’s pressure drop. As 
described above, Forchheimer’s equation must be used to identify 
the characteristics of the porous medium, the viscosity, and inertial 
terms in the two regions. 

This study built an experimental apparatus with a pressure differ¬ 
ence transducer to measure the pressure at the inlet of the canister 


Fig. 5. (a) Experimental image of activated carbon region (without activated carbon 
particles) at a flow rate of 30L/min (open outermost holes on the main sieve 
diaphragm) and (b) velocity vector for simulation of the activated carbon region 
(without activated carbon particles) at a flow rate of 30L/min (open outermost 
holes on main sieve diaphragm). 


outlet ports. Curve fitting is used to obtain Forchheimer’s equations 
for the filter layer and the activated carbon layer. The polynomial 
for the filter layer is 

= auuVs + ft pV s 2 = 27143 Vj + 7849V S 2 . (7) 

L\ 

The polynomial for the activated carbon layer is 

= a 2 (iV s + ftpV 2 = 4274414 + 39284V 2 . (8) 

l 2 

The following coefficients for the filter layer are determined as 
= 1.52 x 10 9 m -2 , ft = 6.41 x 10 3 m _1 . 
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Fig. 7. f-Re' relationships for filter and activated carbon layer. 


Similar coefficients for the activated carbon layer can then be ob¬ 
tained: 

a 2 =2.39 x 10 9 m“ 2 , ft, = 3.21 x 10 4 m _1 . 

The characteristics of the filter layer and the activated carbon 
layer are assumed to be homogeneous and isotropic in CFD simula¬ 
tions. Fig. 6 shows that the simulated pressure drops are very close 
to the experimental values. 

Fig. 7 shows the /-Re/ relationship for the filter and activated 
carbon layers at inhalation rates of 15-135 L/min (inlet Reynolds 
numbers approx. 900-8100). The filter layer exhibits quasi-linear 
(Darcy’s equation) porous medium characteristics. Similarly, the ac¬ 
tivated carbon layer exhibits linear behavior (viscosity effect) at a 
flow rate of 30 L/min. As the inhalation rate increases, non-linear be¬ 


havior (inertial effect) gradually appears. However, the porous flow 
in the activated carbon layer still exhibits a slight inertial effect even 
at a maximum flow rate of 135 L/min. 

3.3. Hole distribution and area of main sieve diaphragm 

For a given thickness of the activated carbon layer, the test models 
were divided into two groups (models Ai, Di, Gi and Hi vs. models 
A 2 , D 2 , G 2 and H 2 , refer to Table 2) to examine the effect of the 
distribution and area of holes in the main sieve diaphragm on the 
pressure drop across the activated carbon layer and flow structure. 
Fig. 8 plots the relationship between the pressure drop and A h /A c 
at a flow rate of 30 L/min (where A h represents the hole area for 
the main sieve diaphragm, and A c represents the cross-section area, 
equal to the area of the main sieve diaphragm). The A h \A c ratio of 
each model is between 0.06 and 0.33, with a value of about 0.09 for 
the prototype canister model Ai and 0.33 for model D 2 . The pressure 
drop for these two groups decreased as the A h \A c value increased, 
revealing that an increase in the hole area can reduce the pressure 
drop for a given hole distribution. The pressure drop for model D 2 
is approximately 113 Pa, which is 30 Pa lower than the 143 Pa value 
for model Ai. The hole areas in models Ai and G 2 are approximately 
equal, but the difference between their pressure drops is 13 Pa. Even 
though the area in model Di with sealed outermost holes exceeds 
that in model A 2 , the pressure drop in model Di is approximately 
5 Pa higher than model A 2 with open outermost holes. This result 
reveals that the effect of the hole distribution on the pressure drop 
of the activated carbon layer exceeds the effect of the hole area. 

Fig. 9(a)-(d) plots velocity contours on an azimuthal angle of 
22.5° plane for models Ai, Di, A 2 and D 2 at a flow rate of 30L/min. 
Since the outermost holes are sealed and the porous flow causes 
momentum loss, the gas flowing through the open holes of the main 
sieve diaphragm does not easily flow past the outer part of the ac¬ 
tivated carbon layer, as Fig. 9(a) and (b) indicates. Thus, a larger 
low-velocity zone forms in this region. Since the hole size of model 
Ai is smaller than that of model Di, the local velocity of gas near 
the open holes in model Ai is also greater than that in model Di. 
The flow structures of these two models are similar except for the 
locations of the open holes on the main sieve diaphragm. Fig. 9(c) 
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Fig. 9. Velocity contours of activated carbon layer of: (a) model Ai, (b) model Di, (c) model A 2 and (d) model D 2 at flow rate of 30L/min. 


and (d) shows that once the outermost holes in models A 2 and D 2 
are opened, the low-velocity region decreases markedly in size. The 
flow structures of the two models are also similar. Fig. 10(a)-(c) 
plots the velocity-position relationship for models Ai, Di, A 2 and D 2 
on a tangential plan view at an azimuthal angle of 22.5°. Fig. 9(a) 
marks three axial positions of the activated carbon layer (z = 0.03, 
0.04 and 0.05 m). Fig. 10(a) shows a non-uniform velocity profile 
for models Ai and D 2 at z = 0.03 m when the outermost holes of 
main sieve diaphragm are sealed. The hole distributions for models 
Ai and D 2 drive the gas passing through the open holes to generate 
a high-velocity zone where the velocity exceeds the mean velocity 
V a ve (^ave = Q/A c ). This forms a low-velocity zone since the gas does 
not flow smoothly into the region near the wall. The hole area ef¬ 
fect leads to a greater local velocity for model Ai than for model Di. 
Models A 2 and D 2 have a more uniform velocity profile than models 
Ai and Di for open outermost holes and a large hole area. The iner¬ 
tial effect of gas passing through the open holes on the main sieve 
diaphragm decreases as the flow distance increases. The local high 
speed flow in models Ai and Di can flow around the outside zone 
of the activated carbon layer at z = 0.04 m (Fig. 10(b)). The velocity 
profiles of these models at the same hole distribution are quite sim¬ 
ilar. Fig. 10(c) reveals that the outlet wall around the corner helps 
form a low-velocity zone near the outlet position (z = 0.05 m) on the 
activated carbon for each test model. All models appear to yield sim¬ 
ilar velocity profiles at the same hole distribution on the main sieve 
diaphragm. 

This study also compares the velocity-position relationships for 
models Ai, Di, A 2 and D 2 at flow rates of 80, 110 and 135L/min 
at z = 0.03 m. A comparison of Fig. ll(a)-(c) with Fig. 10(a) shows 


that the velocity profiles are quite similar in terms of the relative 
positions of the local high and low velocity zones, and the types of 
flow. The flow rate barely affects the velocity profile in the activated 
carbon layer because the porous medium of that layer exhibits quasi- 
linear behavior in this flow range. The high loss of momentum in 
activated carbon layer also causes flow changes in the flow in this 
range that hardly affect the flow behavior, which is dominated by 
the hole distribution on the main sieve diaphragm. 

Figs. 7, 9-11 reveal that the porous flow behavior in the acti¬ 
vated carbon layer does not change markedly at inhalation rates 
from 15 to 135L/min. The local flow structure changes for the test 
models clearly occur in the region near the main sieve diaphragm 
as the gas flows through the holes. The type of passageway on the 
main sieve diaphragm dominates the flow structure in the activated 
carbon layer. Moreover, any given hole distribution yields similar 
flow structures in test models. Forchheimer’s equation for a porous 
medium describes the relationship between the superficial velocity 
and the pressure drop from a macro perspective. A higher super¬ 
ficial velocity corresponds to a greater pressure drop. According to 
the continuity equation, the product of the superficial velocity and 
the flow-through area gives the volumetric flow rate, and the flow¬ 
through area does not equal the cross-sectional area of the activated 
carbon layer. Consider model Ai as an example: the poor hole po¬ 
sitions on the main sieve diaphragm yield a non-uniform velocity 
profile for the intake gas. Gas does not flow easily into the resulting 
low-velocity zone, which creates a dead zone and reduces the flow¬ 
through area. This increases the superficial velocity in the activated 
carbon layer, which in turn increases the pressure drop. In contrast, 
a better hole distribution yields a more uniform velocity profile for 
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a b 



c 




Fig. 10. Velocity profile of models Ai, A 2 , Di and D 2 at 30L/min: (a) z = 0.03m, (b) z = 0.04m and (c) z = 0.05m (0 = 22.5°). 


the intake gas. In model D 2 , for example, increasing the flow-through 
area reduces the size of the dead zone and the superficial velocity 
in the activated carbon layer, thereby reducing the pressure drop. 

In fact, it is difficult to estimate the area of the actual dead zone 
of each model in detail. The pressure drop values in Fig. 8 and 
Eq. (10) can be used to estimate the superficial velocity 14 in the 
activated carbon layer in each model. The flow-through area in the 
activated carbon layer at a constant flow rate of 30 L/min can also 
be calculated. This area is defined as the effective flow-through area 
A e ff. The term A e jy/A c is the ratio of the effective flow-through area 
to the total area, and 1 -A e fflA c is the proportion of the flow-through 
area that is ineffective. Moreover, A c -A e ff represents the area of no 
flow-through, called the dead zone. Fig. 12 shows that model D 2 has 
an effective flow-through area of about 86% and a dead zone of 14%. 
The prototype model Ai has an effective flow-through area of about 
70% and a dead zone of 30%. This wastes the activated carbon layer 
and reveals design flaws. Although the total area of holes in models 
Ai and G 2 is quite similar, the difference in hole distributions re¬ 
sults in different flow structures. The effective flow-through area of 
model G 2 (~75%) is 5% greater than that of model Ai. According to 
the index defined above, the size of the dead zone can be analyzed 
from a macroscopic viewpoint and compared with the CFD simula¬ 
tion results. Using CFD analysis makes it possible to determine the 


effect of the hole distribution and hole area on the flow structure 
and pressure drop. The positions of the dead zones and the prefer¬ 
ential flow can also be identified, which greatly helps improve the 
design of main sieve diaphragms. 


3.4. Activated carbon layer thickness 

This section describes how changing the thickness of the acti¬ 
vated carbon layer affects the pressure drop. Three non-dimensional 
thickness values for the activated carbon layer were simulated-full 
thickness, half thickness, and one-quarter thickness (Table 2). Mod¬ 
els were divided into four groups based on the hole distribution and 
hole area for the main sieve diaphragm. Each group shares the same 
hole distribution and area, but has a different carbon layer thickness. 
Group 1 contains models Ai, Bi and Ci. Group 2 consists of mod¬ 
els Di, Ei and Fi. Group 3 comprises models A 2 , B 2 and C 2 . Group 
4 consists of models D 2 , E 2 and F 2 . Fig. 13 plots the relationship be¬ 
tween the pressure drop per unit length, the superficial velocity, and 
the effective flow-through area ratio of the four groups at an inlet 
flow rate of 30 L/min. Simulation results for all groups show that the 
superficial velocity and pressure drop per unit length increase and 
the effective flow-through area ratio decreases as the carbon layer 
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Fig. 11. Velocity profile of models A u A 2 , Di and D 2 at: (a) 80L/min, (b) llOL/min and (c) 135L/min (z = 0.03 m, 0 = 22.5°). 


170 


160 h 


150 h 


Pm 

<1 


140 h 


130 h 


120 h 


110 


1 1 I 


I 1 1 1 1 I 




I- Groupl (outmost holes closed) 

t- - Group2 (outmost holes opened) 


X c; - 


V 2 


, H , 

X 




I I I 


_L_ 


H o.i 


H 0.095 


H 0.09 


0.085 


H 0.08 


H 0.075 


H 0.07 


0.6 0.65 0.7 0.75 0.8 0.85 

A efl /A 


0.9 



A efl /A c 


0.16 

0.14 

0.12 

0.1 

0.08 

0.06 

0.04 

0.02 

0 


Fig. 12. Relationship between flow efficiency ratio and pressure drop for the eight Fig. 13. Pressure drop per unit length (AP/I) versus superficial velocity (V s ) and 
models. effective flow-through ratio at volumetric flow rates of 30 L/min. 


V s (m/s) 

























C.-C. Li / Chemical Engineering Science 64 (2009) 1832-1843 


1841 


thickness decreases. Within Group 1, in which models Ai, Bi and Ci 
share the same sealed outer holes and original hole area, the super¬ 
ficial velocities greatly increases as the thickness decreases. This re¬ 
sults from the non-uniform velocity profile generated by the sealed 
outermost holes. When the thickness of the activated carbon layer 
decreases, the local high-speed flow will pass quickly through the 
canister outlet. This result causes the gas not flow around the outside 
of the activated and to increase the areas of dead zone. The superfi¬ 
cial velocity of the Group 4 models, with open outermost holes and 
a larger hole area, also increases as the thickness decreases. Their 



Fig. 14. Pressure drop per unit length (AP/I), superficial velocity and effective 
flow-through ratio at volumetric flow rates of 30, 80, 110 and 135L/min. 


more uniform velocity profiles are associated with a lower change 
in the pressure drop per unit length. The difference in the pres¬ 
sure drop per unit length, the superficial velocity, and the effective 
flow-through area for Groups 1, 3 and 2, 4 are attributable to the 
open hole area. For the original hole positions and the hole areas in 
Group 1 models Ai, and Ci at a flow rate of 30L/min, the effec¬ 
tive flow-through area ratios are approximately 0.7, 0.57, and 0.45, 
respectively. Group 4 models D 2 , E 2 , and F 2 have more effective hole 
distributions and larger holes, with effective flow-through area ra¬ 
tios of approximately 0.86, 0.76, and 0.67, respectively. When the 
thickness of the activated carbon layer is decreased to one-quarter 
of its original value, the effective flow-through area and the pressure 
drop per unit length for model F 2 are approximately equal to those 
for model Ai. This reveals that the effect of a more uniform hole 
distribution and a larger hole area on the aerodynamic behavior of 
the flow can be used to reduce the thickness of the activated carbon 
layer. As the thickness of the activated carbon layer is attenuated to 
reduce the distance over which the gas flows, the effect of the hole 
area becomes more apparent. A more suitable distribution of holes 
should be used, and the hole area of the main sieve diaphragm should 
be greatly increased. The effective flow-through areas for test mod¬ 
els at flow rates of 80, 110, and 135L/min are slightly lower than 
those at 30L/min, as Fig. 14 shows. Fig. ll(a)-(c) and Fig. 14 reveal 
that the flow structures of test models, including the velocity pro¬ 
file, the effective flow-through area are barely affected by the flow 
rate, but are significantly affected by the hole positions and area. 

3.5. Air age 

To ensure adequate time for filtration and thus improve the 
canister design, this study considers the residence time of the gas in 
the activated carbon layer together with the pressure drop. The air 
age index was used to investigate aerodynamic behavior within the 
canister. The definition of air age is the local mean age of the fluid, 
or local mean residence time, which is the average time that a fluid 




Fig. 15. Air age contours for activated carbon layer in: (a) model Ai, (b) model Di, (c) model A 2 and model D 2 at Q.= 30L/min, 0 = 22.5°. 
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Fig. 16. Air age distribution for: (a) models Ai, A 2 , Di and D 2 and (b) models Bi, 
B 2 , Ei and E 2 (outlet surface of activated carbon layer, 0 — 22.5°). 


particle takes to reach any point of the domain from the canister 
inlet. The value of air age index is calculated from Eq. (6). The advan¬ 
tage of displacing adsorption time by the air age index is that com¬ 
plex transient adsorption modes must no longer be solved, greatly 
reducing computation time. 

A greater air age within the activated carbon layer represents a 
longer residence time and thus a greater chance for adsorption. How¬ 
ever, a long air age means that the air has difficulty flowing through 
the space, and thus forms a dead zone. Using the flow variable dis¬ 
tributions described by the air age index makes it much easier to 
determine the position of dead zones by observing the flow struc¬ 
ture within the activated carbon layer. The more difficult alternative 
to this approach is to use the velocity profile. 

Fig. 15(a)-(d) shows the air age contours for models Ai, Di, A 2 , 
and D 2 at an azimuthal angle of 22.5° and a flow rate of 30L/min. 


Models Ai and Di, with sealed outermost holes, clearly have a large 
zone of higher air age outside the activated carbon layer (Fig. 15(a) 
and (b)). The gas cannot flow easily through this zone, and a dead 
zone is formed. The position and size of the dead zones in the two 
models are almost identical. Other than the difference between air 
ages near the holes of the main sieve diaphragm, the flow structures 
in other zones for models Ai and Di exhibit similar air-age distribu¬ 
tions since they share the same hole distribution. 

Fig. 15(c) and (d) compares models A 2 and D 2 (open outermost 
holes) to models Ai and Di. The original dead zone apparently 
shrinks after open outermost holes. Models A 2 and D 2 have a larger 
dead zone at the center of the outlet cross-section surface than mod¬ 
els Ai and Di as the fluid flows toward the outer zones. Even if the 
flow structures of models A 2 and D 2 were the same, the different 
hole areas would produce differences in the air age near the holes of 
the main sieve diaphragm. In other zones, the same hole distribution 
creates similar air age distributions. 

Fig. 16(a) plots air age vs. position for models Ai, Di, A 2 and D 2 at 
the outlet surface of the activated carbon layer at an azimuthal angle 
of 22.5°. Model Ai is a commercial product that passes the filtration 
test and has the shortest air age of ~0.67 s. An air age of > 0.67 s 
corresponds to an adsorption time that meets the filtration standard. 
The air ages for models A 2 , Di and D 2 are > 0.67 s and therefore also 
meet the standard. At the outlet of the activated carbon layer, the 
air age distributions for models A 2 and D 2 are greater than those 
for models Ai and Di. The air ages in the middle part are slightly 
larger than ~0.2 s, representing an increase in adsorption time that 
corresponds to a more effective filtration of toxic gases. 

Fig. 16(b) plots air age vs. position for models Bi, B 2 , Ei, and E 2 . 
The air age in most cases is < 0.67 s, indicating that these models 
would fail the standard adsorption time test. Taken together, the 
air age data and the previous computational results identify the to¬ 
tal dead zone and the positions of local dead zones. These finding 
provide a useful reference for designing an optimal main sieve di¬ 
aphragm for gas mask canisters. 

4. Conclusions 

This paper combines flow visualization techniques, pressure drop 
measurements, and CFD tools to investigate a gas mask canister with 
a complex inner configuration of two porous materials. Equations 
for porous media yield parameters that affect the inner flow aero¬ 
dynamics in the canister body, including the distribution and area of 
holes in the main sieve diaphragm and the thickness of the activated 
carbon layer. These findings help pinpoint areas for improvement in 
a prototype canister. 

The simulation results reveal that the channel design of the main 
sieve diaphragm dominates the aerodynamic behavior of the fluid 
within the activated carbon layer. Hole distribution in the main sieve 
diaphragm determines the flow structures of the activated carbon 
layer. Changes in the flow rate do not readily affect the flow struc¬ 
ture. Better hole distribution and a larger hole area correspond to a 
lower pressure drop, a smaller dead zone, and a higher adsorption 
time. However, the structural strength of the main sieve diaphragm 
must help support the weight of the activated carbon layer in highly 
dynamic situations. Deformation that results in a loss of activated 
carbon particles causes the canister to malfunction. Hence, the hole 
position distribution should be optimized by considering not only 
how to reduce the pressure drop and the weight of the main sieve 
diaphragm, but also the structural strength and the limitations of 
the manufacturing technique used. 

The air age can be combined with a defined index, such as the ef¬ 
fective flow-through area ratio (A e jy/A c ), the ineffective flow-through 
area ratio (1 -A e jy/A c ), or the area of the dead zone (. A c -A e ff ), to iden¬ 
tify the size of the entire dead zone, the position of local dead zones 
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and whether or not the adsorption time meets the NIOSH 99.97% 
standard for filtering toxic gases. The results of this study provide 
a valuable reference for designing channels in the main sieve di¬ 
aphragm, and will be helpful in designing gas mask canisters. The pa¬ 
rameter analysis herein reveals that the development of a gas mask 
canister with a low respiratory drop is feasible. 


Notation 

A c cross-section area, m 2 

A e ff flow-through area, m 2 

A h hole area of main sieve diaphragm, m 2 

/ friction factor 

gi gravitational acceleration in the i direction, m s -2 

L porous medium length, m 

p pressure, Pa 

Q. volumetric flow rate, m 3 s -1 

R dimensionless of canister’s radius 

Re, modified Reynolds number 

Reflet inlet Reynolds number 

Si source term which describes the pressure gradient 

caused in porous medium, Pam -1 
U[ velocity component in the i direction, ms -1 

Vave average velocity, ms -1 

V s superficial velocity, ms -1 

Greek letters 


a 

viscous parameter of porous material, m 

P 

inertial parameter of porous material, m 

A p 

pressure drop, Pa 

A z 

air age increase, s 

£ 

turbulent dissipation rte, m 2 s -3 

K 

permeability of porous material, m 2 

Mf 

effective dynamic viscosity, Pa s 

to 

molecular dynamic viscosity, Pa s 

to 

turbulent dynamic viscosity, Pa s 

P 

density, kgnrr 3 

PU'iU'j 

Reynolds stresses, Pa 

or 

turbulent Schmidt number 

T 

air age, s 

T i/ 

viscosity stress tensor, Pa 
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